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We have studied the extended Hubbard model in the atomic limit. The Hamiltonian analyzed 
consists of the effective on-site interaction U and the intersite density-density interactions Wij (both: 
nearest-neighbour and next-nearest-neighbour). The model can be considered as a simple effective 
model of charge ordered insulators. The phase diagrams and thermodynamic properties of this sys- 
tem have been determined within the variational approach, which treats the on-site interaction term 
exactly and the intersite interactions within the mean-field approximation. Our investigation of the 
general case taking into account for the first time the effects of longer-ranged density-density inter- 
action (repulsive and attractive) as well as possible phase separations shows that, depending on the 
values of the interaction parameters and the electron concentration, the system can exhibit not only 
several homogeneous charge ordered (CO) phases, but also various phase separated states (CO-CO 
and CO-nonordered). One finds that the model considered exhibits very interesting multicritical be- 
haviours and features, including among others bicritical, tricritical, critical-end and isolated critical 
points. 

PACS numbers: 71.10.Fd Lattice fermion models (Hubbard model, etc.); 71.45.Lr Charge-density-wave sys- 
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I. INTRODUCTION 

Electron charge orderings phenomena in strongly cor- 
related electron systems are currently under intense in- 
vestigations. Charge orderings (COs) are relevant to 
a broad range of important materials, including mangan- 
ites, cuprates, magnetite, several nickel, vanadium and 
cobalt oxides, heavy fermion systems (e. g. Yb4As3) and 
numerous organic compounds [P-18]. 

Various types of COs have been also observed in 
a great number of experimental systems with local 
electron pairing (for review see [T^HTS] and references 
therein), in particular in the compounds that contain 
cations in two valence states differing by 2e (on-site pair- 
ing) - valence skipping, "negative-?/" centers [17], and 
in the transition metal oxides showing intersite bipo- 
larons e. g. Ti4_a:Va:07, WO3-X, with double charge 
fluctuations on the molecular (rather than atomic) units 
[(Ti4+-Ti''+),(Ti3+-Ti-''+)], etc. 

COs are often found in broad ranges of electron dop- 
ing (e. g. doped manganites ^8nl3j. nickelates [TJ [5], 
Bai_3;Ka;Bi03 |16ffT5] ). In several of these systems many 
experiments showed phase separations involving charge 
orderings [TUHTB] . The CO transitions at T > take 
place either as first order or continuous. Moreover, 
some of CO systems exhibit also a tricritical behaviour 
(e. g. (DI-DCNQI)2Ag [6i). 

An important, conceptually simple model for study- 
ing correlations and for description of charge orderings 
(and various other types of electron orderings) in nar- 
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row energy band systems is the extended Hubbard model 
taking into account both the on-site ([/) and the in- 
tersite {Wij) density-density interactions (the t-U-Wij 
model [roiiraH^). 

In this paper we focus on the atomic limit (tij — 
limit) of the t-U-Wij model. The Hamiltonian considered 
has the following form: 

rrV^ - ~ VFi _ 

— 2^ n,nj - ^ 2^ n„ (1) 

where denotes the creation operator of an electron 
with spin a ='f, J, at the site i, fii = rii„, rii„ — cf^Cia- 
fi is the chemical potential. indicates the 

sum over nearest-neighbour (m = 1) and next-nearest- 
neighbour (to = 2) sites i and j independently. Zm de- 
notes the number of TO-th neighbours. U is the on-site 
density interaction, Wi and W2 are the intersite density- 
density interactions between nearest neighbours (nn) and 
next-nearest neighbours (nnn), respectively. The chemi- 
cal potential /x depends on the concentration of electrons: 

'^=^E('"^»)' (2) 

i 

with < n < 2 and N - the total number of lattice sites. 
{■hi) denotes the average value of the fii operator. 

The model ^ can be considered £lS &I simple model of 
charge ordered insulators. The interactions U and Wij 
can be treated as the effective ones and assumed to in- 
clude all the possible contributions and renormalizations 
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like those coming from the strong electron-phonon cou- 
pling or from the coupling between electrons and other 
electronic subsystems in solid or chemical complexes. In 
such a general case arbitrary values and signs of U and 
Wij are important to consider. 

Notice that the model ([ij can be viewed as the classical 
gas with four possible states at each site and it is equiva- 
lent to a special kind of the Blume-Capel model i. e. the 
S — \ Ising model with single-ion anisotropy and eigen- 
value 'zero' doubly degenerate in an effective magnetic 
field given self-consistently by a value of fixed magneti- 
zation |27l|33ll37|. 

In the analysis we have adopted a variational approach 
(VA) which treats the on-site interaction U exactly and 
the intersite interactions Wij within the mean-field ap- 
proximation (MFA). 

Within the VA the phase diagrams of (ll]) have been 
investigated till now for the case W2 = P5H?7| and the 
stability conditions of states with phase separation have 
not been discussed. Some preliminary results for the case 
W2 7^ have been presented by us in f28]. 

We perform extensive study of the phase diagrams and 
thermodynamic properties of the model ([T|) within VA for 
arbitrary electron concentration n, arbitrary strength of 
the on-site interaction U and the nn repulsion Wi > 0, 
taking into account the effects of interactions between 
nnn W2 (repulsive and attractive). Our comprehensive 
investigation of the general case finds that, depending 
on the values of the interaction parameters and the elec- 
tron concentration, the system can exhibit the charge 
ordered and nonordered homogeneous phases as well as 
(for attractive W2) at least two types of phase separation 
involving charge orderings. Transitions between different 
states and phases can be continuous and discontinuous, 
what implies existence of different critical points on the 
phase diagrams. We present detailed results concerning 
the evolution of phase diagrams as a function of the in- 
teraction parameters and the electron concentration. 

Our studies of the Hamiltonian ([T]) are exact for at- 
tractive W2 in the limit of infinite dimensions. They 
are important as a test and a starting point for a per- 
turbation expansion in powers of the hopping tij and as 
a benchmark for various approximate approaches (like 
dynamical mean field approximation, which is exact the- 
ory for fermion system in the limit of infinite dimensions 
for tij [20|) analyzing the corresponding finite band- 
width models. They can also be useful in a qualitative 
analysis of experimental data for real narrow bandwidth 
materials in which charge orderings phenomena are ob- 
served. 

In the limit W2 — the model (Tj) has been analyzed in 
detail (for review see [27] ^301 and references therein) . In 
particular, the exact solutions were obtained for the one- 
dimensional {d — 1) case (T > 0) employing the method 
based on the equations of motion and Green's function 
formalism (30] or the transfer-matrix method |31| l32] . 
In [33] the phase diagram of as a function of fi for 
W2 = has been derived at T = and confirmed in 133] . 



These studies were based on the Pirogov and Sinai meth- 
ods [3S]. In two dimensional case d = 2 {W2 ^ 0) exact 
ground state diagrams as a function of fi have also been 
obtained [36' using the metod of constructing ground 
state phase diagrams by the reflection positivity prop- 
erty with respect to reflection in lattice planes. 

A number of numerical simulation has also been done 
for W2 — 0. In particular, the critical behaviour near the 
tricritical point have been analyzed using Monte Carlo 
(MC) simulation and MFA [22]. A study of the model in 
finite temperatures using MC simulation has also been 
done for a square lattice (d = 2) [37] [55] . In particu- 
lar, the possibility of phase separation and formation of 
stripes in finite systems was evidenced there. 

In the following we will restrict ourselves to the case of 
repulsive Wi > 0, which favours charge orderings, and 
ziWi > Z2W2- For the sake of simplicity we consider 
mainly two-sublattice orderings on the alternate lattices, 
i. e. the lattices consisting of two interpenetrating sub- 
lattices (every nearest neighbour of every site in one sub- 
lattice is a site in the other sublattice), such as for ex- 
ample simple cubic (SC) or body-center cubic lattices. 
Some preliminary studies of ground state beyond the two- 
sublattice assumption at half-filling (n = 1) are also per- 
formed, which show that for repulsive W2 there is possi- 
bility of occurrence of the multi-sublattice orderings. 

The paper is organized as follows. In section |TT] we 
describe the metod used in this work. There are also 
derived explicit formulas for the free energies of homoge- 
neous phases and states with phase separation as well as 
equations determining the charge-order parameters and 
the chemical potential in homogeneous phases. In sec- 
tion |III| we analyze the properties of the system at zero 
temperature and present ground state diagrams. Sec- 
tion |IV| is devoted to the study of the finite temperature 
phase diagrams for W2 > and W2 < 0. Some partic- 
ular temperature dependencies of the charge-order pa- 
rameter are discussed in section |V] Section I VI I contains 
ground state results for half-filling beyond two-sublattice 
assumption. Finally, section [Vll| reports the most impor- 
tant conclusions and supplementary discussion including 
the validity of the approximation used and the compar- 
ison with real materials. The appendix presents explicit 
expressions of site-dependent self-consistent VA equa- 
tions. 



II. THE METHOD 

The free energy of the system and the self-consistent 
equations for the average number of electrons on sites 
are derived within site-dependent VA in the Appendix. 
Restricting analysis to the two-sublattice orderings the 
explicit formula for the free energy per site obtained in 
the VA has the following form 
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where 



Z„ = l + 2exp[/3(//-Va)] + 
i/jji = nWo + uqWq, tpB 
Wo = ziWi + Z2W2, Wq 



exp[^(/z-2V'„ -[/)], 
= nWo - uqWq, 
= -ziWi + Z2W2, 



and f3 — l/{kBT). The charge-order parameter is defined 
as HQ = (l/2)(n^ - ub), where 11^ = J2iea (^*) ^^^^ 
average electron density in a sublattice a = A, B. 

The condition for the electron concentration Q and 
a minimization of f(n,T) with respect to the charge- 
order parameter lead to a set of self-consistent equations 
(for homogeneous phases): 

(4) 
(5) 



where 



The double 
^ N ^ 



n= (1/2) (n^+vis), 
"Q = (1/2) [nA - nB) , 

{exp [|3i^l - V'a)] + exp [/3(2/i - 2V^„ - U)]} . 



occupancy per site defined 



fii-ffii^) has the following form: 
D = {l/2) (Da + Db), 



as 



(6) 



where Da — exp [/3(2^ — 2ipa — U)]/Za- 

The equations (|4])-(|5]) are solved numerically for T > 
and we obtain nq and fi when n is fixed. The charge- 
ordered (CO) phase is characterized by non-zero value of 
uq, whereas hq = in the non-ordered (NO) phase. 

Let us notice that the free energy ^ is an even func- 
tion of nq so we can restrict ourselves to solutions of 
the set (|4])-(|5]) in the range < nq < 1. It is the result 
of the equivalence of two sublattices. Moreover, (|5| is 
only the necessary condition for an extremum of (|3| thus 
the solutions of (|4])-(|5| can correspond to a minimum or 
a maximum (or a point of inflection) of ([s]). In addition 
the number of minimums can be larger than one, so it is 
very important to find the solution which corresponds to 
the global minimum of (|3|. 

Phase separation (PS) is a state in which two domains 
with different electron concentration: and n_ exist 
in the system (coexistence of two homogeneous phases). 
The free energies of the PS states are calculated from the 
expression: 

fps{n+,n_) = mf+{n+) + (1 - m)f_{n_), (7) 

where f± {n± ) are values of a free energy at n± corre- 
sponding to the lowest energy homogeneous solution and 
m = {n — n_)/(7i_|_ — n_) is a fraction of the system with 
a charge density n+ (n_ < n < n+). The minimization 
of ^ with respect to n+ and n_ yields the equality be- 
tween the chemical potentials in both domains: 



(8) 



(chemical equilibrium) and the following equation (so- 
called Maxwell's construction): 



M+(n+) 
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In the PS states the chemical potential 
^, — ^-^.{nj^) = ^,^{nS) is independent of the elec- 
tron concentration, i.e. d^/dn = 0. 

In the model considered only the following types of PS 
states can occur: (i) PSl is a coexistence of CO and NO 
phases and (ii) PS2 is a coexistence of two CO phases 
with different concentrations. 

In the paper we have used the following conven- 
tion. A second (first) order transition is a transition 
between homogeneous phases with a (dis-)continuous 
change of the order parameter at the transition tem- 
perature. A transition between homogeneous phase and 
PS state is symbolically named as a "third order" tran- 
sition. During this transition a size of one domain in 
the PS state decreases continuously to zero at the tran- 
sition temperature. We have also distinguished a second 
(first) order transition between two PS states, at which 
a (dis-)continuous change of the order parameter in one 
of domains takes place. In the both cases the order pa- 
rameter in the other domain changes continuously. 

Second order transitions are denoted by solid lines on 
phase diagrams, whereas dotted and dashed curves de- 
note first order and "third order" transitions, respectively. 

The phase diagrams obtained are symmetric with re- 
spect to half-filling because of the particle-hole symmetry 
of the hamiltonian ([ij [TSJ [301 EH] , so the diagrams will 
be presented only in the range < n < 1. 



III. THE GROUND STATE 
A. 1^2 > 

In the case of nnn repulsion, i. e. < /c < 1 
{k — Z2W2/Z1W1), the system at T = can exhibit two 
types of CO: ''high CO" (HCO), involving the on-site 
pairing of electrons and "low CO" (LCO), which is the 
ordering without on-site pairs. 

For this case the ground state (GS) diagram de- 
rived within VA as a function of n and U/{—Wq) is 
shown in figure [TJi. At T = HCO can be stable 
phase only if U/i~Wq) < 1 and LCO if J7 > 0. For 
< U/{—Wq) < 1 both types of order can be realized 
depending on n. In the GS, one obtains the following 
results for the ordering parameter ng, the chemical po- 
tential /i and the double occupancy per site D: (i) LCO 
phases: for LCOa (only sublattice A is filled by elec- 
trons without double occupancy ua — 2n, sublattice B is 
empty): nq = n, fj, = 2z2W2n, D = and for LCOb (ev- 
ery site in sublattice A is singly occupied ua = 1, whereas 
ub ~ 2n — I): nq = 1 — n, fi = 2z2W2n — Wq, D = 0, 
(ii) HCO phases (electrons are located only in sublattice 
A, nA — 2n, ub — 0): for HCOa (every site in sublat- 
tice A is doubly occupied): nq = n, fi = U/2 + 2z2W2n, 
D — n/2 and for HCOb (at every site in sublattice A at 
least one electron is located): nq ^ n, ji ~ U + 2z2W2n, 
D = n — 1/2. For n = 1 the transition at U/{—Wq) = 1 
is from the HCO {nq ^ 1, D ^ 0.5) to the NO (Mott 
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FIG. 1. Ground state phase diagrams U/ {—Wq) vs. n: (a) for 
1 > fc > and (b) for A: < {k — —1). Details in text and 
Table U 



TABLE I. PS states occurring in the ground state. 



State 


Domain 


n+ 


Domain 


n_ 


PSIa 


HCO 


1 


NO 









PSIb 


LCCF] 


0.5 


NO 









PSlc 


NO 


1 


LCC 


a 




0.5 


PS2* 


HCO 


1 


LCC 






0.5 



The value of fi in the homogeneous phase is irrelevant here. 



state, HQ — Q, D = Q). In both phases at n = 1 the chem- 
ical potential is ^ U/2 + Wq ^ U/2 + 2Z2W2 - Wq. 
At quarter-fiUing {n = 0.5) in LCO {uq = 0.5, D = 0) 
the value of the chemical potential is /i = U/2 + Z2W2 for 
< U/{~Wq) < 1 and /i = (l/2)Wo for U/{~Wq) > 1. 

For fixed U/{—Wq) the chemical potential /i changes 
discontinuously (except LCO-HCOb and LCO-LCOb), 
whereas for fixed n it is continuous at the phase bound- 
aries. The double occupancy D changes continuously 
at transitions with fixed U/{—Wq), while transitions 
with fixed n are associated with discontinuous change of 
D (except LCO-LCO for n = 1/2 and U/{-Wq) = 1). 
Apart from HCO-NO (for n = 1) and HCOb-LCOb 
transitions, the charge-order parameter uq changes con- 
tinuously in all transitions. 



B. W2<0 

In the GS diagram as a function of n for W2 < (fig- 
ure ^p) one finds also simple linear boundaries between 
various states but in this case the free energies of PS 
states are lower than those of homogeneous phases, apart 
from particular values of concentration: (i) for n = 0.5 
and U/\z2W2\ > 1 the LCO with ng =0.5 is stable, 
(ii) forn = 1 the HCO with ng = 1 if U/{-Wq) < 1 and 
the NO with ng = if U/ (—Wq) > 1 are stable. One can 
also check that the first derivative of chemical potential 
d^/dn in homogeneous phases is negative (apart from the 
ranges mentioned above) what implies that these phases 
are not stable. Definitions of all PS states occurring in 
the GS are collected in Table |T] For U/\z2W2\ <1 only 



the PSIa state (with D = n/2, = U/2 + Z2W2) occurs. 
When 0.5 < n < 1, UI\z2W2\ > 1 and U/{-Wq) < 1 the 
PS2* state with £> = n - 1/2 and ^ = [/ + {3/2)z2W2 is 
stable. For U/\z2W2\ > 1 and n < 0.5 the PSIb with 
D — and fi = Z2W2/2 has the lowest energy, whereas 
for U/\Wq\ > 1 and 0.5 < n < 1 the PSlc with D = 
and fi = {3/2)z2W2 - Wq is stable. 

All transitions in GS for W2 < are associated with 
discontinuous change of the chemical potential fi. The 
double occupancy D changes continuously at transitions 
with fixed U/{—Wq), while transitions with fixed n are 
connected with discontinuous change of D (except LCO- 
LCO for n = 1/2 and [//(-W^q) = 1). 

Notice that for Wi > and W2 < the condition 
U/{-Wq) = 1 implies that U/\z2W2\ > 1, so the line 
U/{-Wq) = 1 is above the line C//IZ2VF2I = 1 and the 
GS phase diagram for any W2 < has always the form 
shown in figure [TJd. 

One should stress that for W2 < the phase stability 
condition is not fulfilled (i. e. dfi/dn < 0) in homoge- 
neous phases except n = 0.5 for U/\z2W2\ > 1 and n = 1. 
It means that the homogeneous phases are not stable. 

For W2 = the free energies of homogeneous phases 
and PS states are degenerated at T = and for such 
a case dfi/dn = in homogeneous phases. This degener- 
ation is removed in any finite temperatures and at T > 
homogeneous phases have the lowest energy. Our dia- 
grams obtained for W2 = are consistent with results 
presented in |27| . 



IV. FINITE TEMPERATURES 

The behaviours of the system for repulsive W2 > and 
attractive W2 < are qualitatively different. 

One obtains from numerical analysis that, for 
ziWi > Z2W2 > taking into account only two- 
sublattice orderings, the PS states are unstable at 
any T > and the homogeneous phases are stable 
{dfi/dn > for any T > 0, even for VF2 = 0). The finite 
temperature phase diagrams have the forms determined 
in [27_, with the replacements: ziWi — >■ —Wq > 0. Thus 
we describe obtained results in this case shortly, direct- 
ing the reader for detailed analyses to [27] . One can con- 
clude that transition temperatures between homogeneous 
phases decrease with increasing W2 > and U/{—Wq). 

For the on-site interaction U/{— Wq) < (2/3) \n2 
and U/{-Wq) > 1 only the second order CO-NO 
transitions occur with increasing temperature. For 
(2/3)ln2 < f//(-VFQ) < 0.62 the first order CO-NO 
transition appears near n = 1 with a tricritical point 
TC connected with a change of transition order 
(e. g. figure [sj;). The TC for 12 — 1 is located at 
ksT/i-WQ) = 1/3 and U / {-Wq) = 2/3ln2. In the 
range 0.62 < U/{-Wq) < 1 the first order CO-CO line 
appears, which is ended at an isolated critical point IC 
of the liquid-gas type (cf. figure |4^). In this case we 
have also a critical end point CE, where three bound- 
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ary lines (one of second order: CO-NO and two of 
first order: CO-NO and CO-CO) connect together. 
The two CO phases: HCO and LCO are distinguish- 
able only in the neighbourhood of the first order line 
CO-CO (HCO-LCO). This line is associated with the 
HCOb-LCOb transition in GS. The first order transi- 
tion CO-NO can exist inside the region 0.79 < n < 1 and 
0.5 < n < 1. The hues consisting of TC, IC and CE 
meet at a new multicritical point, which coordinates are 
n = 0.79, U/{-Wq) = 0.62 and ksT/i-WQ) = 0.24, ap- 
proximately. 



The phase diagrams obtained for attractive W2 < are 
essentially different from those for W2 > 0. The main dif- 
ference is that at sufficiently low temperatures PS states 
are stable. In the ranges of PS states occurrence the ho- 
mogeneous phases can be metastable (if d^i/dn > 0) or 
unstable (if d^/dn < 0). In the homogeneous phases oc- 
curring at higher temperatures (above the regions of PS 
occurrence) the stability condition dji/dn > is fulfilled. 



For clarity of the presentation we will discuss the be- 
haviour of the system at T > for W2 < distinguish- 
ing three regimes of U/{—Wq): on-site attraction (sec- 
IV A I , strong on-site repulsion (section |IV B ) and 



tion 



weak on-site repulsion (section IV C\ 



We also distinguish three different PSl states and four 
different PS2 states (labeled by subscripts A, B, C or 
superscript *). All states in a particular group are states 
with the same type of phase separation (i. e. CO-NO or 
CO-CO), however they occur in different regions of the 
phase diagram and such distinction has been introduced 
to clarify the presented diagrams (cf. especially figure |5|. 
A similar distinction has been done for all critical points 
connected with the phase separation (four S-type points: 
B, B', B", B*, three T-type points: T, T', T" and several 
points of H-, E-, F- types). The lines M-N-O and X-Y 
indicate the first order transitions between PS states on 
the kBT/{—WQ) vs. n diagrams (cf. figures [s] and |4| . 



A. The case of on-site attraction 

For any on-site attraction {U < 0) the phase diagrams 
are qualitatively similar, all (second order and "third or- 
der") transition temperatures decrease with increasing U 
and for U = the transition temperatures account for 
a half of these in the limit U —00. In a case of W2 < 
beyond half-filling the PS states can be stable also at 
T > 0. Examples of the ksT / {—Wq) vs. n phase dia- 
grams evaluated for U/{—Wq) = and various ratios of 
k = Z2W2/Z1W1 < are shown in figure[2) If < |fc| < 1 
the (homogeneous) CO and NO phases are separated by 
the second order transition line. 

When -0.6 < fc < (figure [2^) a "third order" transi- 
tion takes place at low temperatures, leading first to the 
PS into two coexisting CO phases (PS2a), while at still 
lower temperatures CO and NO phases coexist (PSIa)- 
The critical point for this phase separation (denoted as 
B, we shall call this point a bicritical endpoint, BEP) 
is located inside the CO phase. The E-F solid line (we 
shall refer to i?-point as a critical endpoint, CEP) is asso- 
ciated with continuous transition between two different 
PS states (PS1a-PS2a, the second order CO-NO tran- 
sition occurs in the domain with lower concentration). 

For k < —0.6 (figure [2]d) the transition between PS 
states does not occur, the area of PS2a stability vanishes 
and the critical point for the phase separation (denoted 
as T, which is a tricritical point, TCP) lies on the sec- 
ond order line CO-NO. As k —>■ —00 the T-point occurs 
at n = 1 and the homogeneous CO phase does not exist 
beyond half-filling. 

When k = —0.6 the lower branch of the "third or- 
der" curve approaches the critical point (H) paraboli- 
cally. The tricritical behaviour for k < —0.6 changes into 
the bicritical behaviour for k > —0.6. The i/-point is 
a higher order critical point (HCP) and in this point the 
lines consisting of B E, F, and T points connect to- 
gether (for fixed U/{—Wq)). Similar scenario takes place 
also for the on-site repulsion apart from the bicritical be- 
haviour connected with the PS2* state, which can exist 
also for k > —0.6). 
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n n 

FIG. 3. (Color online) Phase diagrams kBT/(—WQ) vs. n for Wi > and several values of U/{~Wq) and k = Z2W2/Z1W1: 
(a) U/{-Wq) = 0.2, k = -0.2; (b) U/{-Wq) = 0.2, k = -0.4; (c) U/{-Wq) = 0.6, k = -0.2; and (d) U/{-Wq) = 0.6, k = -1. 
Dotted, solid and dashed lines indicate first order, second order and "third order" boundaries, respectively. 



One should notice that a type of the critical point for 
separation (which can be BEP, TCP or HCP) is modi- 
fied only by a change of the strength of the nnn attrac- 
tion. The location of the transition lines between ho- 
mogeneous phases on the A;bT/(— Wq) vs. n diagrams is 
not affected by the value of W2 (the transitions are at 
the same ksT / (—Wq) as in the previous case of W2 > 0, 
which only depends on the on-site interaction U for fixed 
n) . Effectively, the transition temperatures increase with 
increasing strength of attractive W2 ■ 

The labels of the critical points for phase separations 
are given with a correspondence to those in |29j . where 
the Ising model with nn and nnn interactions was con- 
sidered. Notice that our model is equivalent to the Ising 
one in the U — ^ ±cx) limits. 



B. The case of strong on-site repulsion 

For any U/{—Wq) > 1 the structure of the phase di- 
agrams and the sequences of transitions are similar as 
those in the previous case (for corresponding values of 
fc), but now the double occupancy of sites is strongly 
reduced due to repulsive U and the phase diagrams are 
(almost) symmetric with respect to n = 0.5 (cf. figureflb 
and table|ll. B', H', T', E' and F' points (as well as , 
H", T", W and F" points) appear (cf. |28]), which corre- 
spond to B, H, T, E and F points, respectively. Critical 
behaviours at A' and A" points are the same as at A 



points {A^ B,H,T). 

The exact symmetry occurs at [/ — > +00. In this 
limit the phase diagrams are the diagrams for U — > —00 
with re-scaled axes: /^^^/(-VFq) -> /cb7'/(-4M^q) and 
n -> (l/2)n. 

Transition temperatures are only weakly dependent on 
the on-site repulsion and one finds a small decrease of 
them with increasing U for U/{—Wq) > 1. 

C. The case of small on-site repulsion 

The range < U/{—Wq) < 1 is the most interesting 
one and the phase diagrams are more complicated than 
those in previous cases. Due to the variety of the be- 
haviour in this regime of on-site interaction we only 
present some particular examples of the phase diagrams 
(figures |3] and |4]). 

In this range of U/{—Wq) an occurrence of the PS2* 
state for 0.5 < n_ < n < n+ < 1 is possible and the crit- 
ical point for this phase separation (BEP, denoted as 
B*) is located inside the CO phase also for |fc| > 0.6. 
It contrasts with BEPs mentioned previously (i. e. B, 
B' and B"), which occur only for |fc| < 0.6. The re- 
gion of PS2* occurrence extends from the ground state if 

U/\Z2W2\ > 1. 

When < U/{-Wq) < 2/3 In 2 one can see on the 
phase diagrams a new behaviour, i. e. discontinuous tran- 
sitions between PS states connected with the M-N-O 
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FIG. 4. (Color online) Phase diagrams fcsr/(— Wq) vs. n for Wi > and several values of U/{~Wq) and k = Z2W2/Z1W1: 
(a) (7/(-VKq) = 0.8, fc = -0.2; (b) U/{-Wq) = 0.8, fc = -0.8; (c) (//(-VFq) = 0.94, k = -0.4; and (d) U/{~Wq) = 0.94, 
fc = — 1. Dotted, solid and dashed lines indicate first order, second order and "third order" boundaries, respectively. Near 
n — 1, on the right of the PS2* occurrence region, the CO phase is stable, what is not shown on the diagrams (c) and (d) 
explicitly. 




FIG. 5. (Color online) Phase diagrams fc_Br/(— Wq) vs. U/{~Wq) for several values of n and k = Z2W2I z\W\: (a) k — —0.2, 
n — 0.25; (b) k — —0.2, n = 0.75; (c) k = —1, n — 0.25; and (d) k — —1, n — 0.25. Dotted, solid and dashed lines indicate first 
order, second order and "third order" boundaries, respectively. One should notice that some axes are broken. 
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line. On figure [SJd two such transitions are presented, 
i. e. PSIa-PSIb and PS1a-PS2*. There are also pos- 
sible PS2a-PS2b and PS2a-PS2* transitions, which are 
not shown there. All transitions between homogeneous 
phases are second order in this range. 

In the range 2/31n2 < U/{-Wq) < 1, besides the be- 
haviours mentioned previously, the critical points con- 
nected with transitions between homogeneous phases 
can appear on the phase diagrams: (i) TC for 
U/{-Wq) < 0.62 and the first order CO-NO boundary 
near n = 1 (cf. figures [sj: andjsji) and (ii) CE and IC for 
U/{-Wq) > 0.62 with the first order HCO-LCO tran- 
sition slightly dependent on the electron concentration 
(cf. figures |4]) . It leads to an appearance of isolated ar- 
eas of PSlc state stability (only for k < —0.6, figure |4|3) 
and to discontinuous transitions between two PS states: 
(i) PS2*-PSlc (figures ^ and ^) and (ii) PS2*-PS2c 
(only for < 0.6, cf. also figure ^jp) in a restricted 
ranges of U/{—Wq). These behaviours are connected 
with the X-Y line, which is linked with the first order 
border lines CO-NO and CO-CO. For > fc > -0.6 
above the X-Y line a second order PSlc-PS2c transi- 
tion also occurs. The temperature associated with X-Y 
line is independent of k and depends only on U/{—Wq ). 

The finite temperature phase diagrams as a func- 
tion of U/{—Wq) at fixed n for k = —0.2 and k = —1 
are shown in figure [s] The first order boundaries: 
(i) PSIa-PSIb and PS2a-PS2b (on the diagrams for 
n = 0.25) and (ii) PS1a-PS2* and PS2a-PS2* (on the 
diagrams for n = 0.75) are associated with the M-N-O 
line, whereas the first order boundaries: (iii) PS2*-PSlc 
and PS2*-PS2c (on the diagrams for n — 0.75) are con- 
nected with the X-Y lines. For k — —0.2 and n = 0.25 
at higher temperatures the PSIa and the PSIb are not 
distinguishable and on the diagram (figure [5j:) the first 
order boundary line ends at a critical point of the liquid- 
gas type (similar to IC). In figure [sjl for n = 0.75 one can 
also see following sequence of transitions with increasing 
temperature: PS2*-C0-PSlc-C0-N0. It is interest- 
ing to notice that the PSlc state exists here at higher 
temperatures than the homogeneous CO phase (see also 
figure |4j3) . 

One should notice that the "third order" boundaries 
in figure [5] are not the Hues of BEPs (nor TCPs). For 
considered ranges of U /{—Wq) and k the n-coordinates 
of these points fulfill the following conditions: (i) T and B 
points: n < 0.5, (ii) T and B' points: n < 0.25, (iii) T", 
B" and B* points: n > 0.75. One should also remember 
that E-F, E'-F', E"-F", M-N-O and X-Y lines in their 
ranges of occurrence are independent of n. 



V. CHARGE-ORDER PARAMETER VS. 
TEMPERATURE 

In this section we present two representative tempera- 
ture dependencies of the charge-order parameter for fixed 
model parameters. 



(a) '■" 
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FIG. 6. (Color online) Temperature dependencies of the 
charge-order parameter uq for (a) U/{—Wq) = 0.2, n = 0.2 
and k = -0.4; and (b) U/i-Wg) = 0.94, n = 0.75 and 
k = -0.4. 



In figure [6^ we have plotted the charge-order parame- 
ter UQ as a function of ksT/i-WQ) for U/{-Wq) = 0.2, 
n = 0.2 and k = -0.4. At A:bT/(-VFq) 0.058 one ob- 
serves a discontinuous change of uq > in one domain 
(in the other uq = 0), connected with the first-order 
transition PSIa-PSIb- For kBT/i-Wg) = 0.073 the 
continuous transition PS1b-PS2b occurs and uq raises 
to non-zero value in the domain with lower electron con- 
centration (cf. figure |3]3). At kBT/i-Wq) = 0.079 the 
domain with higher electron concentration vanishes con- 
tinuously, a "third order" transition PS2b-C0 occurs and 
the whole system is characterized by one value of the 
charge order-parameter. Next, at kBT/{—WQ) = 0.167 
Uq goes to zero, i. e. one has the second order CO-NO 
transition. 

Finally, let us comment on the temperature depen- 
dence of ng for U/{-Wq) = 0.94, n = 0.75 and k = -0.4 
(figure At kBT/i-Wg) = 0.043 one observes 

a first order PS2*-PSlc transition, which is connected 
with discontinuous change of uq in the domain with 
higher electron concentration (cf. figure |4|;). Next, at 
kBT/i-Wq) = 0.069 a continuous PSlc-PS2c transi- 
tion occurs (now one has continuous change of ng in the 
domain with higher electron concentration). At higher 
temperatures one can notice a "third order" PS2c-C0 
transition (at A;bT/(— Wq) = 0.075) and a second order 
CO-NO transition (at kBT/i-Wg) = 0.196). 



VI. BEYOND THE TWO-SUBLATTICE 
ORDERINGS FOR REPULSIVE W2 

The nn repulsion Wi > 0, as well as the nnn attrac- 
tion W2 < 0, favour two-sublattice ordering and in such 
a case no other types of long-range order can occur. 
On the other hand, the repulsive W2 > compete with 
Wi > reducing stability of two-sublattice orderings and 
can yield the appearance of multi-sublattice orderings. In 
this section we consider charge-orderings on the regular 
lattices taking into account not only two-sublattice or- 
derings. We will consider them as follows. 

More general, we can define the charge-order param- 
eter as riq- = exp (ig • i?i/a), where determine 
a location of i site in the space and a is a hypercubic 
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lattice constant. 

The four-sublattice orderings can be considered on the 
alternate lattices, in which both interpenetrating sublat- 
tices are also alternate lattices. Examples of such lattices 
are ID-chain and 2D-square (SQ) lattice. 

When we consider four-sublattice orderings in the GS, 
the homogeneous phases not mentioned in section [llT| are 
found to occur for W2 > in some definite ranges of k 
and U/ziWi. 

For example at half-filling a so-called island charge or- 
dered phase (ICO, . . . 2200. . . , in d = 1) or stripe charge 
ordered phase (SCO, in d = 2, q— (0,7r)) can occur for 
Z2W2 / ziWi > 0.5. In HCO phase we have g = tt in d = 1 
and q — (tt , tt) ill d = 2. The GS phase diagram taking 
into account the four-sublattice orderings on 2D-square 
lattice for n = 1, Wi > is shown in figure |7^. In the 
case of ID-chain the range of ICO occurrence is the same 
as that of SCO in 2D. In both cases the VA results (for 
n = 1) are consistent with exact results |331IMlf5SlHD . 

On the SC lattice we cannot consider four-sublattice 
orderings, because the two interpenetrating sublattices 
are fee sublattices and they are not alternate lattices. In 
case of such a lattice at half-filling the following three 
types of commensurate charge orderings should be con- 
sidered: (i) q = (0, 0, tt) (plane charge ordered phase, 
PCO), (ii) g=(0,7r,7r) (SCO) and (iii) g=(7r,7r,7r) 
(HCO). The GS phase diagram for SC lattice and n = 1 
is shown in figure [7J3. The PCO phase does not occur in 
GS and a region of the NO phase stability is extended in 
comparison to the lower dimension cases. 

In all CO phases mentioned previously the number of 
electrons on the particular site can be = or rii = 2 
and the charge order parameter in each phase has a max- 
imum possible value, i. e 
(rij = 1 at every site). 

One should notice that a region of the HCO phase oc- 
currence (on the k vs. U/{ziWi) diagram) does not de- 
pend on the lattice dimension. This result is in agreement 
with GS phase diagrams obtained in section |III[ The 
discontinuous transition HCO-NO is at A; = 1 — U/ziWi 
what is equivalent to U/{~Wq) — 1 (for k < 0.5 and 
71 = 1). 

Let us stress that we have not analyzed the four- 
sublattice orderings at T > 0. They can be stable at 
sufficiently low temperatures, such as SCO near half- 
filling for k > 0.5 lU or some states with phase sepa- 
ration between different CO phases or between CO and 
NO phases (for n ^ 1) |42]. Thus, the finite temperature 
phase diagrams for W2 > with taking into account the 
four-sublattice orderings can be in general more involved 
than those discussed in section II VI 



n^= 1. 



In the NO Ug- = 



VII. CONCLUDING REMARKS 

In this paper we studied atomic limit of the ex- 
tended Hubbard model with intersite nn repulsion Wi. 
By taking into account for the first time the effects of 




(b) 



SCO 

(0,7t,7t)/ 



NO 
(0,0,0) 



HCO 



FIG. 7. Ground state phase diagrams for n — 1 with consid- 
eration of the multi-sublattice orderings: (a) for 2D-square 
lattice (for ID-chain the region of SGO is replaced by that of 
ICO), (b) for SC lattice in d = 3. The phases are also labeled 
by vector g (fc = 22W2/21VF1, Wi > 0). 



nnn density-density interaction (attractive and repulsive 
W2 ^ 0) and including into consideration the states with 
phase separation (involving CO) our paper substantially 
extends and generalizes the results of previous works con- 
cerning the model considered |26| [57] . Let us summarize 
the most important conclusions of our work: 

(i) Depending on the values of the interaction param- 
eters and the electron concentration, the system can ex- 
hibits not only several charge ordered and nonordered 
homogeneous phases, but also (for attractive W2) var- 
ious phase separated states involving charge orderings: 
PSl(CO-NO) and PS2(C0-C0). 

(ii) Obtained phase diagrams have a very rich struc- 
ture with multicritical behaviours. The regions of PS 
states (both PSl and PS2) stability expand with increas- 
ing of the next-nearest-neighbour attraction. Moreover, 
the transitions between PS states can be both continuous 
and discontinuous such as those between homogeneous 
phases. 

(iii) The value of W2 < determines a type of the 
multicritical point associated with PS states (which we 
mentioned as bicritical, tricritical or high- order critical 
point). 

(iv) For repulsive W2 > there is a possibility of oc- 
currence of multi-sublattice orderings, e. g. stripe (or 
island) structures. In particular, these types of ordering 
have the lowest energy at T = for n = 1 (cf. section VI I 
if fc = Z2W2/Z1W1 > 0.5 and U/ziWi > 1 as well as for 
n = 0.5 if fc > 0.5 and U/ziWi — > -l-oo and in such cases 
they can be stable also at sufficiently low temperatures. 

Let us stress that for Wi > and W2 < the derived 
results are exact in the limit of infinite dimensions, where 
the MFA treatment of intersite interactions becomes the 
rigorous one. 

Charge ordered phases are stable (for W2 — 0) when 
the quantum perturbation by finite bandwidth is intro- 
duced mmSHlS]. Thus one can conclude that the PSl 
and PS2 states, which involve the CO phases, should also 
occur in the presence of hopping term. The stability of 
the phase separated CO-NO state in finite temperature 
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(for tij and W2 = 0) in a definite range of the electron 
concentration was confirmed by using dynamical mean 
field approximation [5U] . If tij ^ for [/ < the on-site 
superconducting states can occur, whereas for J7 > it is 
necessary to consider also magnetic orderings |16l HBH^ . 
In such a case various phase separation states involving 
superconducting and (or) magnetic orderings can also be 
stable. 

For the case Wi < 0, which was not analysed in the 
present work, the model ([ij {W2 — 0) exhibits a phase 
separation NO-NO (electron droplet states) at low tem- 
peratures |50l I51| . In this PS state different spatial non- 
ordered regions have different average electron concen- 
trations. In such a case, at higher temperatures only the 
homogeneous NO phase occurs. 

The VA results can be exact in the limit of infinite 
dimensions only. Below we discuss briefiy some results 
obtained for W2 = within other approaches. In par- 
ticular, the Bethe-Peierls- Weiss (BPW) treatment of 
the Wi term predicts the following ranges of the exis- 
tence of long-range CO at T = for SQ and SC lat- 
tices in the limits: (a) U/ziWi — >■ -|-cxi (LCO phase): 
l/zi<n< (2zi - l)/zi and (b) U/ziWi -00 (HCO 
phase): 2/zi < n < 2{zi — l)/zi. In particular: 

(i) d = 2 (SQ lattice): (a) 0.25 < n < 0.75 and 
1.25 <n< 1.75 (LCO), (b) 0.5 < n < 1.5 (HCO); 

(ii) d = 'S (SC lattice): (a) 0.17 < n < 0.83 and 
1.17 <n< 1.83 (LCO), (b) 0.33 <n< 1.67 (HCO); 

(iii) d = +00 (hypercubic lattice): (a) < n < 1 and 
1 < n < 2 (LCO), (b) < n < 2 (HCO). 

The Monte Carlo calculations performed for SQ lat- 
tice [37; yields for the case U +00 even more re- 
stricted ranges for long-range CO (LCO) at T = 0: 
0.37 < n < 0.63 and 1.37 < n < 1.63. In this particular 
case the existence of percolations of the effective clusters 
has been confirmed. These percolations vanish at transi- 
tion temperature. 

The phase diagrams for Wi > and W2 = obtained 
by exact solution for the Bethe lattice |51] (which 
is equivalent to BPW approximation) have a similar 
structure as the VA diagrams (even for small zi = 3). 
The main difference is a reentrant behaviour found in 
the case of Bethe lattice for U <0 if 71 < 2/zi and 
n > 2{zi — l)/zi, where the sequence of phase transi- 
tions: NO->CO^-NO can occur with increasing tempera- 
ture. The transition temperatures determined in [51' are 
in general lower than those obtained in VA, but obviously 
in the limit of large coordination number the rigorous re- 
sults for Bethe lattice reduce to those of VA. 

Comparing the GS diagram obtained for W2 = in 
VA (figure [iji) with the exact one for ID-chain [3^ , we 
notice that all the border lines are the same, although 
in the exact solution for d = 1 the long-range charge or- 
derings in GS exist only in the ranges (i) 0.5 < 71 < 1 
and < U/{-Wq) < 1, (ii) n = 0.5 and U/{-Wq) > 0, 
(iii) n = 1 and U/{—Wq) < 1 (what corresponds to the 
regions of HCOb, LCO and HCO phases existence in fig- 
ure flk, respectively) . The values of fi and D obtained 



in VA are consistent with exact ones for arbitrary n and 
U/{—Wq). Moreover, the GS phase diagrams as a func- 
tion of /i for W2 = derived within VA agree exactly 
with the corresponding rigorous solutions in d = 1 and 
d = 2 I33H3S]- 

The above discussion implies that VA in the case 
W2 = can give qualitatively reasonable results beyond 
the percolation thresholds also for lattices of finite di- 
mensionality and this statement should also be true for 
W2 ^ 0, at least for small attractive W2. 

The electron concentration n and chemical potential 
fi are (thermodynamically) conjugated variables in the 
bulk systems fSO]. However, one can fit the concentra- 
tion rather than the potential in a controlled way ex- 
perimentally. In such a case /i is a dependent internal 
parameter, which is determined by the temperature, the 
value of n, and other model parameters (cf. (|2|). Thus 
the obtained phase diagrams as a function of the con- 
centration are quite important because in real systems 
n can vary in a large range and charge orderings are of- 
ten found in extended ranges of electron doping (e. g. in 
doped manganites [8,-13, , nickelates [TJ |5] and the doped 
barium bismuthates [TBl - fT5] ). In Bechgaard salts the con- 
centration is n = 1/2 [SHI]. In charge transfer salts n 
changes, dependent on the pressure, in the vicinity of 
n — 2/3, whereas for several complex TCNQ salts n is 
near n = 1/2 [S5]. In cuprates [T^ and refer- 
ences therein) and in conducting polymers [53] n is near 
half-filling in the insulating state and it strongly changes 
under doping. 

Although our model is (in many aspects) oversimpli- 
fied, it can be useful in qualitative analysis of experi- 
mental data for real narrow-band materials and it can 
be used to understand better properties of several CO 
systems mentioned above and in section [T] 

In particular, our results predict existence of the phase 
separation (CO-NO, CO-CO) generated by the effective 
nnn attractive interactions and describe their possible 
evolutions and phase transitions with increasing T and 
a change of n. The electron phase separation involving 
COs is shown experimentally in several systems quoted 
above, e. g. in Ri„2;Ca2;Mn03 (R=La, Bi, Nd, etc.), at 
dopings ranging from x = 0.33 to x = 0.82 ^MIEl EQ]. 
Among the materials for which the on-site local electron 
pairing (valence skipping) has been either established or 
suggested (cf. section |l| the best candidates to exhibit 
the phase separation phenomena are the doped barium 
bismuthates (BaPbi-jjBi^rOa and Bai-^jKaBiOa) |16ffT51 
HH] . For these systems, being oxide perovskites, a very 
large dielectric constant strongly weakens the long-range 
Coulomb repulsion, which is the main factor preventing 
the phase separation [ST. 

Our results show that also the transitions at T > be- 
tween various homogeneous phases (HCO and LCO) and 
nonordered states can be either first order or continuous 
ones and both these types of the CO transitions are ex- 
perimentally observed in real narrow-band materials [lOf - 
il3j . Moreover, the theory predicts that with a change of 
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the model parameters {U/{—Wq), W2, n) the system can 
exhibits various types of multicritical behaviour (includ- 
ing TCP, BCP, etc.) resulting from the competition of 
the on-site repulsion {U > 0) and the effective intersite 
repulsion Wq < 0. In fact, some of charge ordered sys- 
tems are found to exhibit the multicritical behaviour, e. g. 
in organic conductor (DI-DCNQI)2Ag (T^ = 210 K) the 
temperature vs. pressure phase diagram shows continu- 
ous and first order boundaries with a tricritical point [B] . 
The increasing pressure changes first the order of transi- 
tion, resulting in a tricritical point, then it yields a com- 
plete suppression of charge orderings at any T. 



A is defined as (A) = 



Tri? means a trace 



of any operator B and it is calculated in the Fock space. 

The explicit formula for the free energy obtained in the 
VA has the following form 



1 



In Zj 



(A.3) 



where 



Z, = l + 2exp [Pin - V,)] + exp - 2^;, - t/)], 
while the expression for the average number of electrons 
at z-site is given by 
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Appendix: Site-dependent self-consistent VA 
equations 

Within the VA the on-site interaction term is treated 
exactly and the intersite interactions are decoupled 
within the MFA (site-dependent): 

h.fi] {fii) fij + {iij) hi - {hi) (hj) . (A.l) 

A variational Hamiltonian for the model ([ij has a form 

1 



Uhi-fhii ~ fi.fii - -^n-ii^i 



(A.2) 



where ipi = J2j^i ^'.j'^ji A^i = M - "0* and rii = {hi). Hq 
is diagonal in representation of occupancy numbers (i. e. 
Hi is diagonal in the base consisting of |0), | t)i I i)- 
I ti) at i site with eigenvalues: 0, doubly degenerated 
—Hi, and U — 2fii, respectively) and a general expression 
for the free energy F in the grand canonical ensemble in 
the VA is 
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Tr 



exp(-/3i/o)] } + f^{Ne) 



where /3 — Ne — J2i (^e) — the number 

of electrons in the system. The average value of operator 



n, = —{exp - + exp [|3{2^l - 2^-, - U)]}, 

(A.4) 

so one has a set of iV -I- 1 self-consistent equations to solve 
consisting of N equations in form of ( |A.4[ ) (for every site 
from iV sites) and the condition (pi in the form: 



= 0. 



(A.5) 



The double occupancy Di of the site i is determined 
by the following equation: 

A = {n^tn^i) - ^ exp [/3(2^ - 2i^, - U)] . (A.6) 



The solutions of the set ( |A.4 )-( A.5 1 can correspond to 
a minimum, a maximum or a point of inflection of the 
free energy (A.3) on the (iV — 1) -dimensional manifold 
in A^-dimensional space \ni\f^i defined by the condition 
(A.5 1. To find the solutions corresponding to stable (or 



metastable) states of the system, one should find a con- 
ditional minimum of F with respect to all rii with the 



condition (A.5 1 



One can prove that for two-sublattice orderings, if 
F = F{nA, ns) have a conditional minimum with respect 
to and with condition (A.5 1, then F = Fining) 
have also a minimum with respect to uq = (jia — nB)/2 
(if n is fixed). So the procedure used in section does 
not lose stable (metastable) solutions. In this instance, 
F is the free energy only of homogeneous phases (as 
an assumption) and cannot describe any phase separated 
states, which energies are calculated from ([t]). One needs 
to check the stability condition djjL/dn > for the ho- 
mogeneous phases, which is one of the sufficient con- 
ditions for the conditional minimum of F. In the case 
of two-sublattice orderings on the alternate lattices the 
Eqs. (|A.3| -( |A.4[ ) reduce to Eqs. obtained in sec- 

tion [nj 
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